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We report the observation of the baryonic B decay B° — > A% AK~ with a significance larger 
than 7 standard deviations based on 471 x 10 6 BB pairs collected with the BABAR detector at the 
PEP-II storage ring at SLAC. We measure the branching fraction for the decay B° —>■ At AK~ to 
be (3.8 ± 0.8 s tat ± 0.2 sys ± 1.0, + ) x 10~ 5 . The uncertainties are statistical, systematic, and due to 
the uncertainty in the A£ branching fraction. We find that the A+K~ invariant mass distribution 
shows an enhancement above 3.5GeV/c 2 . 

PACS numbers: 13.25.Hw, 13.60.Rj, 14.20.Lq 



While baryons are produced in (6.8 ± 0.6)% [l[ of all 
B-meson decays, little is known about the detailed me- 
chanics of these decays and more generally about hadron 
fragmentation into baryons. We can increase our under- 
standing of baryon production in B decays by comparing 
decay rates for related exclusive final states. In this paper 
we present a measurement of the decay B° — > A+AK~ 
0. Currently, no experimental results are available for 
this decay. 

This analysis is based on a dataset of about 429 fb -1 , 
corresponding to 471 x 10 6 BB pairs, collected with the 
BABAR detector at the PEP-II asymmetric-energy e + e~ 
storage ring, operated at a center-of-mass energy equal 
to the T(AS) mass. The signal efficiency is determined 
with a Monte Carlo simulation based on EvtGen Q for 
the event generation, and GEANT4 [4] for the detector 
simulation. The B° — > A+AK~ Monte Carlo events are 



generated uniformly in the A+AK~ phase space. Monte 
Carlo simulated events are used to study background con- 
tributions as well. 

The BABAR detector is described in detail elsewhere 
Q . Charged particle trajectories are measured by a five- 
layer double-sided silicon vertex tracker and a 40-layer 
drift chamber, both immersed in a 1.5 T axial magnetic 
field. Charged particle identification is provided by ion- 
ization energy measurements along with Cherenkov ra- 
diation detection by an internally reflecting ring-imaging 
detector (DIRC). 

The A+ is reconstructed in the decay mode A£ — > 
pK~TT + and the A in the decay mode A — > pn + . For 
the identification of proton, kaon, and pion candidates, 
we use selection criteria based on the measurements of 
the specific ionization in the tracking detectors, and of 
the Cherenkov radiation in the DIRC IfJ. 
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For the identification of the p coming from the 
the average efficiency is about 95% while the probabil- 
ity of misidentifying a kaon as a proton is less than 2%. 
The average efficiency for the K ~ identification is about 
90% The probability of misidentifying a pion as a kaon 
is about 5%. These are the dominant misidentification 
probabilities for each particle type. The A£ daughters 
and the A daughters are each fit to a common vertex and 
the A+ and the A candidate invariant mass is required 
to lie within 3<r of the world average mass[H; i.e., in the 
range 2.273 to 2.299 GeV/c 2 and 1.113 to 1.119 GeV/c 2 , 
respectively. For the reconstruction of the B candidate, 
the mass of the A£ candidate is constrained to its nom- 
inal value [l[ and is combined with a A and a K~ can- 
didate. Since the A candidate mass is already well mea- 
sured, it is not constrained. 

The A+, A and K~ candidates are then fitted to a com- 
mon vertex and the confidence level of this fit is required 
to exceed 0.2%. 

A possible source for fake signal events is the decay 
B° — > A+pK~ir + which has the same final state as 
the decay under investigation. In order to suppress this 
background we require that the distance between the B 
vertex and the A vertex in the xy plane (with z par- 
allel to the beam axis) exceeds 0.4 cm. This constraint 
reduces combinatoric background by 18%, and the back- 
ground from B° —> A+pK~n + by 99.6%. The expected 
remaining background from this decay is determined to 
be 0.1 ± 0.1 eventsQ. 

The separation of signal and background in the 
candidate sample is obtained by using two kine- 
matic variables, AE — E* B — s/s/2 and toes = 

\J (s/2 + Pi • p B ) 2 /Ef - \p\ 2 B , where y/a is the e + e~ 
center-of-mass energy and E* B the energy of the B can- 
didate in the center-of-mass system. {Ei, pi) is the four- 
momentum vector of the e + e~ system and ps the B- 
candidatc momentum vector, both measured in the lab- 
oratory frame. For true B decays toes is centered at the 
B-meson mass and AE is centered at zero. B candidates 
are required to have an to-es value between 5.272 and 
5.288 GeV/c 2 . 

Figure[T]shows the AE distribution of the selected can- 
didates, fitted in the range from —0.12 to 0.30 GeV. We 
fit the signal with a Gaussian with the mean \i and width 
cr fixed to the values obtained from a fit to the Monte 
Carlo simulation (p = 0.247MeV and a = 8.381 MeV), 
leaving only the signal yield floating. The background is 
described by a first-order polynomial. A binned maxi- 
mum likelihood fit with this probability density function 
(PDF) gives a signal yield of 51 ± 9 events. (For the 
branching fraction measurement described later, we use 
the excess number of candidates above background as 
the estimate of the number of signal events.) The con- 
fidence level for the null hypothesis, considering statisti- 
cal uncertainties only, is 2.6 x 10~ 15 , which corresponds 



to a statistical significance of 8 standard deviations. A 
possible background from B° —> A+XJ°K~, which rises 
slowly up to AE sa —0.06 GeV and drops sharply be- 
tween AE = -0.05 GeV and -0.02 GeV, is not visible in 

Fig.m 
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FIG. 1: The AE distribution for A^AK~ candidates in data 
with all selection criteria applied (points). The solid line is 
the overall fit result, while the dashed line is the background 
component. 

Since the decay dynamics of baryonic B decays are 
largely unknown, we investigate the invariant-mass dis- 
tribution of the two-body systems. Intermediate states 
would appear as differences in the invariant-mass distri- 
bution for date and for the B° -> A+AK~ Monte Carlo 
simulation, in which the final state is generated accord- 
ing to three-body phase space. Using the same function 
that we used in Fig. [TJ we fit the AE distributions for 
ten ranges of the three two-body masses. The results are 
compared to the phase space Monte Carlo simulation in 
Fig. [2] While the m(A+A) and m(AK~) distributions 
show no significant deviations, the m(A+K~) distribu- 
tion shows the data concentrated in the upper half of the 
allowed mass range, contrary to the Monte Carlo simula- 
tion. A possible explanation for this is a resonant decay 
via a baryon resonance that has not yet been observed. 
Another possibilty is enhanced rates at both m(A^A) 
and m(AK~) thresholds. 

Because the efficiency varies over the Dalitz plot and 
the distribution of candidates in data is unknown a pri- 
ori, we must use the distribution of the data events in the 
Dalitz plot to estimate the efficiency. The small number 
of candidates, combined with resolution and edge effects, 
make the simple weighting of events by the inverse of 
the efficiency problematic. Instead we determine a set 
of weights to apply to the simulated events so that the 
resulting weighted Monte Carlo distributions mimic the 
data. We make the assumption that the dependence of 
the decay dynamics on the two-body invariant masses 
can be factorized into the product of three functions that 
each depend on one invariant mass and weight the Monte 
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TABLE I: Two-body invariant mass distributions used for the 
weighting and the corresponding reconstruction efficiency e. 
The dots indicate how often the respective mass distributions 
are used to determine the weights. 

m(A+K~) m(A+A) m(AK~) Efficiency e in % 

10.90 
8.60 

• 9.21 

• • 9.19 

• • 8.81 
•• • 8.80 
•• •• 8.77 
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FIG. 2: The m(yl+yl), m(AK') and m(AtK~) distributions 
in data (points) in comparison with the Monte Carlo sample 
(red histogram). 



Carlo events with the function: 

w[m(A+K~), m(A+A), m(AK~] 

= w a [m(A+K-)} ■ w b [m{A+A)] ■ w c [m{AK-)}. (1) 

By dividing the background subtracted m(AfK~) dis- 
tribution (Fig. [2J by the corresponding distribution 
from the phase-space Monte Carlo simulation, we ob- 
tain the weights w a [m(A+ K ~)], which are used to weight 



the Monte Carlo candidates. (If a negative weight 
is required the weight is constrained to zero.) Next, 
we use the weighted Monte Carlo candidates to de- 
termine Wb[m(A+A)} in the same way and then use 
w a {m(A+K~)] ■ Wb{m(A+A)] to determine w c [m(AK ~)]. 
After each weighting we determine the reconstruction 
efficiency by a fit to the AE distribution for weighted 
Monte Carlo. Starting with these weights, the weight- 
ing is repeated until the reconstruction efficiency con- 
verges and the two-body mass distributions in data and 
Monte Carlo agree within statistical uncertainties. The 
efficiency after each weighting is shown in Table |TJ Since 
the m(A+K~) distribution in data shows the strongest 
deviations compared to the phase space Monte Carlo sim- 
ulation we use the efficiency e obtained after the second 
weighting in m(A+K~) [e = 8.81%]. The comparison 
between data and weighted Monte Carlo events in the 
two-body masses can be seen in Fig. [3] Note that, by con- 
struction, the data and simulation agree exactly for the 
m(A+K~) distribution. The close agreement in the other 
two distributions shows that the form given in Eq. (|TJ) is 
adequate to describe any correlations between variables 
in the data. The effect of the statistical uncertainties 
in the data on the efficiency determination is described 
below. 

For the branching fraction calculation we determine 
the number of reconstructed events by a fit to AE with 
a first-order polynomial for the background. To avoid a 
potential bias introduced by an assumption of the signal 
shape, we fit the region -0.12 < AE < 0.30 GeV exclu- 
sive of the signal region -0.03 < AE < 0.03 GeV. By 
extrapolating the background yield into the signal region 
and subtracting it from the integral of the histogram in 
this region we obtain a signal yield of N s i g = 50 ± 11. 
This results in a branching fraction of 
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FIG. 3: The background-subtracted m(A^A), m(AK~) and 
m(A£ K~) distributions for data (points) and for the final 
weighted Monte Carlo sample (red histogram). 



with N B g = N So +N B o = (471±3) x 10 6 , assuming equal 
production of B°B° and B+B~ in the decay of the Y(4S). 
The branching fractions B(A+ -> pK-ir+) = (5.0±1.3)% 
and B(A ->B7r + ) — (63.9 ± 0.5)% are the world averages 
from Ref. 

Several sources of systematic uncertainties are inves- 
tigated and summarized in Table |TTJ Most of the un- 
certainties are derived from comparisons between Monte 



Carlo simulations and control samples in data. System- 
atic uncertainties arise from uncertainties in charged par- 
ticle reconstruction efficiencies (0.9%) and charged par- 
ticle identification efficiencies (2.4%), and from statisti- 
cal uncertainties in the Monte Carlo simulation (0.5%). 
The systematic uncertainty on the number of BB pairs is 
0.6%. The systematic uncertainty from the A branching 
fraction amounts to 0.8%. 

The systematic uncertainty introduced by neglecting a 
possible B° — > A+E°K~ background is determined by 
adding a PDF for this background to the fit function 
used for the AE fit shown in Fig. [TJ Allowing nonnega- 
tive contributions from this background, the fit returns 
a value of 0.0^0 o- For a conservative limit on this sys- 
tematic uncertainty we fix the yield to 1.8 and take the 
change in the number of signal events as systematic un- 
certainty (1.0%). The AE distribution in Fig. [1] shows 
an enhancement below — 0.14GeV, caused by decays of 
the type B — > A+AK~ir. Due to the limited resolution 
these events could leak into the fit region from —0.12 to 
0.30 GeV. We determine the resulting systematic uncer- 
tainty by changing the fit region to —0.10 to 0.30 GeV. 
The branching fraction changes by 1.8%. 

The uncertainty arising from the chosen background 
description is determined by repeating the fit to deter- 
mine the signal yield with a second-order polynomial for 
the background. The number of signal events changes 
by 2.0%. A comparison between data and Monte Carlo 
events shows that the mean of the AE distribution in 
data is shifted by —0.003 GeV. We determine the re- 
sulting systematic uncertainty by shifting the signal re- 
gion for the fit in the Monte Carlo AE distribution by 
0.003 GeV. This changes the efficiency to 8.60%, corre- 
sponding to a systematic uncertainty of 0.8%. 

The uncertainty due to the treatment of the three-body 
phase space in the efficiency correction is estimated from 
the variation of the efficiency when performing further 
iterations of weighting. Table U shows a variation from 
8.81% down to 8.77%. This corresponds to a systematic 
uncertainty of 0.5%. 

A possible additional contribution to the statistical 
uncertainty coming from the efficiency determination, 
where we determined the weights based on data events in 
ranges of the two-body masses, is studied by performing 
the efficiency correction in ranges of m(A+K~) only, with 
unweighted Monte Carlo events. The change in over- 
all reconstruction efficiency is negligible compared to the 
statistical uncertainty. 

Adding all contributions in quadrature we obtain a 
systematic uncertainty of 4.1%. The significance of the 
signal, including additive systematic uncertainties is de- 
termined to be more than 7 standard deviations. This 
significance includes systematic uncertainties from the 
B° A+S°K- and B A+AK~tt background as well 
as the AE background description. 

In summary, we observe the decay B° — > A^AK~ with 
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TABLE II: Summary of the relative systeamtic uncertainties 
on the branching fraction B(B° — ¥ A+AK~). 



Source 


Relative syst. uncert. 


Track reconstruction 


U.y/o 


Charged particle ID 




N BB 


0.6% 


B(A -> P7T+) 


0.8% 


Monte Carlo statistics 


0.5% 


B° A+E°K~ 


1.0% 


B -> A+AK-% 


1.8% 


AE background description 


2.0% 


fj,{AE) shift 


0.8% 


Efficiency determination 


0.5% 


Total 


4.1% 



a significance larger than 7 standard deviations and mea- 
sure a branching fraction of 

B(B° -> yl+lK") 

= (3.8 ± 0.8 stat ± 0.2 sys ± 1.0^+) x 10- 5 . (3) 

The decay rate is not uniform over three-body phase 
space; rather, it is dominant at high A~^K~ mass. 
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